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Abstract 

Cardassian models are an alternative to general relativity which have 
been proposed as an approach to explaining accelerated cosmic expansion 
while avoiding directly introducing dark energy. They are generally for- 
mulated only in the homogeneous and isotropic case. In this paper an 
extension of the usual formulation to general spatially homogeneous ge- 
ometries is given. A characteristic feature of many classes of Cardassian 
models is the occurrence of big rip singularities where the scale factor 
tends to infinity after a finite time. It is shown that big rip singulari- 
ties are also widespread in more general homogeneous cases. It is also 
shown that there is isotropization in the approach to a big rip singularity 
which bears a strong resemblance to the late-time isotropization observed 
in cosmological spacetimes which accelerate forever in the future. 



1 Introduction 

The realization, based on observational data, that the expansion of the universe 
is accelerated tias led to a lot of work which modifies traditional cosmological 
models by introducing exotic types of matter or replacing general relativity 
by something else. The focus of interest in the following is an example of the 
second approach. It concerns the Cardassian models which were first introduced 
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in [3]. The standard solutions of the Einstein equations apphed in cosmology 
are those which are homogeneous and isotropic. These are the Friedmann- 
Lemaitre-Robertson- Walker (FLRW) models. Often attention is concentrated 
on the spatially flat case, where the spacetime metric takes the form 

- dt^ + (t) {dx^ +dy^ + dz^) (1) 

with scale factor a{t). In these solutions a central role is played by the Friedmann 
equation 

^\"^^JLP + ^ (2) 



where p is the energy density, A is the cosmological constant and geometrized 
units are used where the gravitational constant and the speed of light take the 
numerical value unity. The other independent Einstein equation is the evolution 
equation 

a 47r , , , A 

ad 6 
where p is the pressure. As a consequence of energy-momentum conservation 
the matter quantities p and p satisfy 

p + 3H{p + p) = (4) 

where H = d/a is the Hubble parameter. It is well known that if equation ^ is 
satisfied at some time to and equations ([3]) and ([4]) are satisfied at all times then 
([2]) is satisfied at all times. This is a special case of the general phenomenon of 
propagation of the constraints in solutions of the Einstein equations. Acceler- 
ated expansion means that a > and it follows from ([3|) that this is only possible 
if A > or p -I- 3p < 0. Many matter models in general relativity satisfy the 
strong energy condition p + 3p > which corresponds to the fact that gravity is 
attractive. Accelerated expansion is only possible in general relativity if either 
A > or there is a matter field (dark energy) which violates the strong energy 
condition. The idea of [3] is to replace general relativity by another theory so 
as to make accelerated cosmological expansion possible with a vanishing cos- 
mological constant and without the need for matter violating the strong energy 
condition. 

In the Cardassian models of [3], which are formulated for homogeneous and 
isotropic geometries, the Friedmann equation ([2]) is replaced by the equation 

d\ ^ Stt 



, Ap + fip)) (5) 

Q> / o 

for a smooth function /. The standard choice for f{p) has the power-law form 
Bp" for constants B and n but in principle the formalism can be developed for 
any function /. In later work 4J this approach was changed and the Cardassian 
term is only a function of part of the energy density of the matter, the rest 
energy density. This splitting requires the input of more information than the 
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energy-momentum tensor alone. The analysis in this paper will concentrate on 
the original version although the later version will be commented on. 

In [3] a possible origin of the additional term in the Cardassian framework 
was discussed in the context of braneworld models. A Cardassian term also 
arises in loop quantum cosmology as a manifestation of the discreteness of space- 
time [TT]. In the latter case the constant B is negative. 

If the Cardassian models are to be regarded as defining a theory of grav- 
ity then there should exist a definition of these models applying to solutions 
without symmetry which reduces to the known definition in the homogeneous 
and isotropic case. A general theory of this type is also of practical importance 
for comparisons with cosmological observations. Some of these comparisons use 
linearization around the homogeneous and isotropic models. In order to ob- 
tain the field equations for linearized perturbations the natural procedure is to 
start with a general theory and linearize. To the authors' knowledge no such 
general theory has been found. The comparison of theory with some aspects 
of observations, e.g. estimates of the age of the universe or supernova data, 
only requires the scale factor but other data such as that on cosmic microwave 
background fluctuations requires more structure (including an understanding of 
linear perturbations) on the theoretical side to allow a comparison. The fact 
that calculations could be done for linearized perturbations is based on the use 
of what will be called the 'effective fluid' in what follows. 

Because of the high symmetry of homogeneous and isotropic geometries any 
Cardassian model with a specific choice of / and matter modelled by a perfect 
fiuid with a specific choice of equation of state p — g{p) is equivalent to general 
relativity with an effective fluid with equation of state p — g{p), where g is 
determined uniquely by / and g. Starting from this observation it is possible to 
use the linearized Einstein-Euler equations with this effective equation of state 
to do perturbation theory. If, on the other hand, another specific description of 
matter, such as kinetic theory, is used then things are more complicated. For a 
particular solution of the Einstein equations coupled to this matter model which 
is homogeneous and isotropic it is possible to define energy density and pressure 
and, with a suitable monotonicity assumption, an equation of state relating the 
two. However this 'equation of state' depends on the solution (and hence on the 
initial data) and as a consequence not only on the choice of matter model. Thus 
it cannot be used to define a generalization to the inhomogeneous case. The 
unfortunate lack of a general theory cannot simply be overcome by an appeal to 
braneworld models or loop quantum cosmology since these constructions have 
also only been carried out fully under symmetry assumptions. 

It will be shown that it is possible to define a natural generalization of the 
Cardassian models to spacetimes which are homogeneous but not necessarily 
isotropic. This accomodates any form of matter. For general matter models 
the energy-momentum tensor is intrinsically anisotropic - it has distinct prin- 
cipal pressures - and so this cannot be replaced by any effective perfect fluid 
description. A limitation of this deflnition is that, as will be shown in section 3, 
for the most general classes of homogeneous geometries and matter it requires 
the assumption that the cosmological expansion never vanishes in order to have 
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a regular evolution equation. An interesting dynamical feature of Cardassian 
models in contrast to general relativity with conventional fluid is the occur- 
rence of big rip singularities This means that the scale factor a{t) tends 
to infinity in a finite proper time. The main result of this paper is that under 
general assumptions homogeneous spacetimes isotropize in the approach to a 
big rip singularity. This is analogous to the cosmic no hair theorem for the late- 
time behaviour of solutions of the Einstein equations with positive cosmological 
constant and its proof in the homogeneous case due to Wald , 13J . 

The structure of the paper is as follows. Section 2 contains an analysis of 
the dynamics of homogeneous and isotropic Cardassian models using an effective 
potential. The method of an effective fluid is also described, together with its 
limitations. In section 3 the generalization to cases which are homogeneous 
but not isotropic is presented and isotropization near a big rip singularity is 
proved. In the last section possible further developments of Cardassian models 
and related topics are discussed. This paper is based in part on the diploma 
thesis of the first author [T]. 



2 Isotropic models 

The basic equation for Cardassian models is The analysis of this section 
is restricted to isotropic models but spatial curvature is included so that (O 
becomes 

a\ Stt , _ k , . 

a) + 

The curvature parameter /c is -)- 1 , or — 1 . The equation ([4]) expressing energy- 
momentum conservation remains unchanged. Differentiating ([6]) with repect to 
time and using Q in the result gives the following evolution equation: 



47r , . Stt 

-y(. + 3,) + - 



f(.p)-lip + p)f'ip) 



(7) 



It is clear from this expression that accelerated expansion is possible even when 
the matter satisfies the strong energy condition. These equations for describing 
gravity can in principle be combined with any model of matter. The equations 
of motion of the matter are assumed to take exactly the same form as in general 
relativity. The remainder of this section will concentrate on the case where the 
matter is described by a perfect fluid with equation of state p — g{p). Other 
matter models will be included in the discussion of anisotropic geometries in 
the next section. 

The dynamics of isotropic models with a fluid will now be analysed. A 
convenient way of doing this is to use the analogy of the motion of a particle in 
a potential. It is assumed that the function g defining the equation of state is 
non-negative with g{0) = and satisfies < g'{p) < 1 for all values of p. Define 
a quantity r representing the mass density of the fluid by 

r{p)^exp\ [\a + g{a))-'da\ (8) 
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Under the given assumptions < dr/dp < 1. From ^ it can be concluded that 
r < p < for p > 1 and < p < r for p < 1. It follows that the mapping 
from p to r can be inverted. Integrating equation ^ shows that a^r is time 
independent. Combining this with the conservation law results in a relation 
p = h{a) for a function h depending on g. In the particular case that g is a 
linear function of p the function h is equal to a constant times a power of a. 
Define a potential y by F (a) = —^a'^{h{a) + f{h{a))). Then the dynamics of a 
is equivalent to that of a particle in this potential. The particle has energy — fc/2 
as a consequence of The picture with the potential is useful for giving an 
intuitive method for seeing how the dynamics looks. It can however do more and 
give rigorous proofs of the main features of the dynamics. For the convenience 
of the reader this will be shown in detail in an appendix. 

Consider the case, which may be called the standard case, where the Car- 
dassian term is defined by f{p) — Bp" and the equation of state of the perfect 
fluid is g{p) = (7 — l)p with 1 < 7 < 2. Then p ~ Aa^^^ for a positive constant 
A and 

V{a) = -_(^a-37+2 + BA"a-3»7+2)^ (9) 
3 

The cases n = and n — 1 will be excluded from consideration since the first can 

be identified with general relativity with cosmological constant and the second is 

general relativity with vanishing cosmological constant in disguise. These cases 

are already sufficiently well known. The domain of definition of the potential is 

(0, 00) since A is positive and —87 + 2 is negative. The derivative V'{a) has a 

or _B < and n > 
37 37 



zero a, if and only if cither B > and n < or B < and n > In both 



BA"-^{2 - 3n7) 



(10) 



37-2 

Important information about the potential is provided by the sign of V" at the 
critical point. 

y"(a,) = -^a,-37^[(_3^^2)(-37+l)+BA"-i(-3n7+2)(-3n7+l)a;='^""'^''] 

(11) 

From this it follows that for n > 1 the potential has a minimum at a* while for 
n < 1 it has a maximum there. It is now easy to read off information about the 
qualitative behaviour of the solution using the results of the Appendix. When 
doing this Figs. [T]-[2]are useful for the purposes of orientation. In each of these 
figures a representative potential is shown for the given range of parameters. 
Since the universe is expanding today only those solutions will be considered 
which arc expanding at some time. Note that this excludes the exceptional 
time-independent solutions which occur when B > 0, n < 01 B < 0, n > 
and k — V{a^,). The results will be summarized in a proposition. 

Proposition 1 Let a(t) be the scale factor of an isotropic and homogeneous 
solution of a Cardassian model with the matter being a perfect fluid with linear 
equation of state p ~ (7 — and 1 < 7 < 2. Suppose that there is a time to 
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Figure 1: The cases B > 0,n< and B < 0,-^ < n< 1. 



Figure 2: The case B < 0, n > 1. 



with d{to) > 0. Then if neither B nor n — ^ is zero the following mutually 
exclusive cases occur: 



1. The solution is periodic; occurs for all solutions with B < Q, n > 1 and 
k>0. 

2. a starts from zero at a finite time, increases to a maximum and returns to 
zero in finite time; occurs for all solutions with B > Q, n < F(a*) > 

—k/2 and a < a^:, all solutions with B > 0, n > and k > 0, all solutions 
with B <0, j^<n<l,F(a*)>— fc/2 and a < and all solutions with 

B < and n < ^ 



3. a tends to oo in the distant past, is decreasing before a certain time, reaches 
a minimum value and is increasing thereafter; occurs for all solutions with 
B > 0, n < —k/2 < V{a*) and a > a*, all solutions with B < Q, 
< n < 1, ~k/2 < y(a*) and a > a^, and all solutions with B < 0, 

n > 1 and fc < 

4- a increases in a monotone way from zero at a finite time; all solutions 
with B > 0, n < and —k/2 > ^(a*), all solutions with B > 0, < ?^ 
and k <0 and all solutions with S<0, ^<n<l and —k/2 > y(a*). // 
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— fc/2 > V{a^) then a oo as t ~> oo while if — fc/2 = V{a^) the solution 
approaches a* as t oo. 

The case 5 = corresponds to the weU-known case of general relativity. 
If n = ^ the solutions behave as in the second or fourth items in the list of 

the proposition, depending on the value of i^i^i^ in comparison to —k/2. 

They are equivalent to models without Cardassian term and with a modified 
value of the curvature. One of the most interesting types of behaviour included 
in the statement of the proposition is the fact that for certain cases the big 
bang singularity is replaced by a bounce. It is this which, in particular, makes 
oscillating solutions possible. 

Next more detail will be obtained about the leading order asymptotics of 
the scale factor as a approaches zero or infinity. Without loss of generality it 
may be assumed that the approach of a to zero takes place in the past and its 
approach to infinity takes place in the future. Consider first the limit a —> 0. Let 
the degenerate case when n — I and i? = — 1 be excluded from consideration. 
For a small V{a) = — iC^a~^P(l + o(l)) as a — > for constants C and p with 
p > 1/2. When a — > the curvature parameter is negligible compared to the 
potential and so a = Ca~P{l + o{l)). Moreover one of the terms in the potential 
dominates the other. Hence 

a = {C{p+l))^{t-ta)^il + o{l)) (12) 

where to is the time at which a vanishes. If n > 1 the Cardassian term dominates 
as a — > and putting in the relevant constants gives 

a= {12nBn^j'^A'')^{t-to)^il + o{l)). (13) 

If n < 1 the Cardassian term is negligible for a small and 

{12'Kj'^A)^{t-to)^{l + o{l)). (14) 

At a big rip singularity, where a tends to infinity at a finite time i» it is 
again the case that the curvature parameter becomes negligible and one of the 
terms in the potential dominates the other. In this case a = Ca'?(l + o(l)) 
with g > and C > 0. For a finite-time singularity it is necessary that q > 1. 
This corresponds to n < 0. Letting b = a^^ and reversing the direction of time 
recovers the situation analysed above, with a weaker restriction on the powers. 
In fact q — 1 corresponds to p = 1. Putting in the relevant constants allows the 
following asymptotics to be derived using the previous analysis: 

a= il27rBn^j^A'')^{t^ - t)^ (1 + o(l)). (15) 

Here t^ is the time at which the singularity occurs. With some more work 
the o(l) error terms in these expressions can be replaced with specific positive 
powers oi t — to OT t^, — t repectively [l] . 

Now the effective fiuid approach to isotropic Cardassian models will be de- 
scribed. The idea is to define a new energy density by p = p + f{p) and a new 
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pressure p so that dp/dt = —^{p + p). If this can be done and if the expression 
for p in terms of p can be inverted a relation of the form p = g{p) can be ob- 
tained. Then the dynamical properties of the Cardassian model are equivalent 
to that of the Einstein equations coupled to a perfect fluid with equation of 
state p = g{p). It should, however, be noted that the effective fluid has exotic 
properties such as violating the strong and weak energy conditions in general, 
the latter in the case B < 0. In this construction the relation 

P = P + fip)ip + p)-f{p) (16) 

is obtained. To solve for p in terms of p it suffices for the condition f'{p) > — 1 
to be satisfied since then the inverse function theorem can be applied. In the 
case of a power-law Cardassian term this works in the case Bn > but not when 
Bn < 0. In the latter case the attempt to carry out this construction leads to an 
equation of state which is not single-valued. Thus this transformation cannot 
be carried out globally for Bn < 0. Note that in the case of forever expanding 
models for which p — > as i ^ oo the transformation can be carried out locally 
so as to apply to solutions at late times. When f{p) = Sp" and p = {'y — l)p the 
equation of state for the effective fluid approaches the linear form p = (717 — l)p 
as p — ^ 0. By comparison with models with a fluid with exactly linear equation 
of state it may be guessed that power-law inflation will occur for n in the range 
(^, 1). This is conflrmed by the formula (fT3|) . 

To end this section the alternative approach to Cardassian models will be 
considered briefly. In this case the expression p + f{p) in ^ is replaced by 
Pi + f{pi) + P2 for two independent quantities pi and p2- Here pi represents the 
rest energy density and p2 all other forms of energy contained in the matter. 
The equations of motion for the matter are dpi/dt — ~3Hpi and dp2/dt = 
—SH(p2 + p)- In the case that p = (7 — l)p the dynamics can be modelled by 
the potential 

A-jr 

Via) = (^la^i + BA>~3n+2) ^ A2a~^^+^. (17) 

3 

The critical points of this potential are not so easy to determine as those of 
the potential in ([9]) and it may not be possible to compute them explicitly. 
Nevertheless it should be possible to do a complete analysis of the qualitative 
features of the dynamics with some more work. 

3 Isotropization near big rip singularities 

There is a natural way to generalize the Cardassian models to the case of space- 
times which are homogeneous but not isotropic. The basic principle is to re- 
quire the propagation of the constraint equations. The Hamiltonian constraint, 
which corresponds to the Friedmann equation ^ in the isotropic case, is mod- 
ified when passing from general relativity to the Cardassian case in exactly the 
same way as was done for isotropic models - the energy density p is replaced 
by p -f /(p) and all other terms are left unchanged. The momentum constraint 
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is left as it is. The equations of motion of the matter, and hence the equation 
for energy-momentum conservation, are as in general relativity. Then the ques- 
tion of the propagation of the constraints is posed. If both the constraints are 
satisfied at some time t = tg and a certain evolution equation generalizing ([7]) 
is satisfied at all times together with the equations of motion of the matter are 
the constraints satisfied at all times? It turns out that there is an evolution 
equation which gives a positive answer to this question and reduces to the usual 
evolution equation for general relativity when the Cardassian term / vanishes. 
This equation will now be written down. 

The spatially homogeneous spacctimes can be divided into the Bianchi space- 
times of types I-IX and the Kantowski-Sachs spacetimes. In the following only 
Bianchi models will be considered, although it would presumably be possible 
to carry out a similar analysis in the Kantowski-Sachs case. In the case of the 
Bianchi models it is possible to assume without loss of generality that the spatial 
manifold is simply connected. For passing from a Bianchi model to its univer- 
sal cover does not change the dynamics. Once this has been done the spatial 
manifold can be identified with a simply connected three-dimensional Lie group 
G and the metric and other geometrical objects can usefully be parametrized 
by their components in a basis of left-invariant one-forms 9'' on G with a dual 
basis Ci of vector fields. The spatial metric is specified by a matrix gij{t) de- 
pending only on time. The spacetime metric is of the form —dt'^ + QijO^O^ . Let 
kij be the second fundamental form of the hypersurfaces of constant t and de- 
fine H = — ifrfc. In the isotropic case this quantity H reduces to the Hubble 
parameter introduced previously. 

The modified Hamiltonian constraint is 

R - %fc'^ + (trfc)2 = 16^(p + f{p)). (18) 

The evolution equation for the second fundamental form is 

dk"^ ■ 1 

= R'j + (trfc)Fj - 87r(S'*j - -trSS]) - inpS} + Wj (19) 

with trace 

rUirlc\ 

= R+ {tvkf + AntrS - Unp + trU. (20) 



dt 

Combining this with the generalized Hamiltonian constraint gives 

= i(trfc)2 + 47r(p + tr5) + trU + 16^/(p). (21) 
Here Sij = Tij is the spatial projection of the energy-momentum tensor and 



-2fip) + f'{p) (^p 



tik trk 



5). (22) 



In the isotropic case ((2T|) and ((22)) reduce to ([7]). Note that j is in general 
only well defined as long as trfc ^ since there are two terms with tik in the 
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denominator. The first of these vanishes in Bianchi Class A while the second is 
regular in the case of an untilted perfect fluid. (For the definitions of these terms 
the reader is referred to [H].) Without these restrictions it could happen that 
the evolution breaks down. It is however the case that for the types of solutions 
which will be of most interest in the following some control of the apparently 
singular terms can be obtained. The source of this control are some inequalities 
implied by the Hamiltonian constraint which will now be presented. Equation 
can be rewritten as 

^(trfc)2 = a,y^ ~R+ 16^(p + f{p)) (23) 

where aij is the tracefree part of kij. In Bianchi types I- VIII it is known that 
R < [13]. Thus if a solution is of one of these Bianchi types and if / is non- 
negative then if trfc tends to zero at some time it follows that p tends to zero, 
with p < 23^(tr/c)2, f{p) < j^itrk)^ and a.^a'^ < |(trfc)2. 

In spacetimes which are not isotropic the momentum constraint becomes 
non-trivial and must be taken into account. As has already been mentioned, 
for the Cardassian models it is assumed that the momentum constraint takes 
exactly the same form as in general relativity. This assumption is used in the 
calculation to verify that the Hamiltonian constraint propagates. In showing 
that the momentum constraint itself propagates the key feature of the modified 
evolution equation which is used is that the modification only affects the trace 
while the tracefree part of the equation remains unchanged. 

It will now be investigated how solutions of these equations behave near a big 
rip singularity. This means by definition that the determinant det g of the spatial 
metric or equivalently, in more geometrical language, its volume form tends to 
infinity in finite time. Assumptions will be made which are analogous to those 
made in [T3] . Consider a spacetime of Bianchi type I- VIII which is expanding 
at some time and satisfies the Cardassian equations introduced above. Suppose 
further that the matter content of spacetime satisfies the dominant and strong 
energy conditions. Thus there is no dark energy included. Finally, suppose that 
the Cardassian term is of the form f{p) = Bp" with B > and n < 0. 

It will now be shown that in a spacetime of this type trfc remains negative 
as long as a smooth solution exists. In fact it will be shown that as long as trfc 
is bounded there cannot exist a sequence of times i„ approaching the end of 
an interval where the solution exists for which trfc(t„) — > 0. For this purpose, 
suppose that there is a smooth solution defined on a time interval [to^T) with 
trfc bounded and a sequence t„ —>■ T with the property in question. Let Y 
denote the right hand side of (pij) . Then 

trfc(i) = trfc(i„) - Y{s)ds (24) 

Using the remarks following ([23]) it can be seen that all terms in Y except tr[/ 
can be bounded by an expression of the form C(trfc)^ for a constant C. Under 
the assumption that trfc is bounded this implies that trU < C|trfc| for a constant 
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C. The terms in trU involving / satisiy the same type of bound. It remains to 
consider the expressions in Y where trfc occurs in the denominator. In order to 
obtain an estimate for Vij* it is useful to evaluate it in an orthonormal frame. 
There the components of j can be bounded by p due to the dominant energy 
condition while the rotation coefficients are of order one. Thus |Vij*| < Cp for 
a constant C and |Vij'|/(trfc) < C|trfc|. To deal with the last term in iiU first 
note that 

fc',„/(trfc) = a'„/(trfc) + ^5'^. (25) 

Thus there is an explicit contribution ^f'{p)trS to trt/. This is no more difficult 
to estimate than the terms already discussed. For the remaining quantity to be 
estimated it is possible to use the inequality that a^jSh < {S^jSlY^'^^a^ jO^ i)^/"^ . 
The dominant energy condition shows that any component of S in an orthonor- 
mal frame is bounded by p. Hence the last expression can be bounded by Cp. 
Putting all this information together and writing t — tn —t gives an inequality 
of the form ^ 

|trA:|(r) < |tr/c(0)| + C / \trk{s)\ds. (26) 
Jo 

Applying Gronwall's inequality [5] and letting n —^ oo gives the conclusion that 
trfc = on the whole interval, in contradiction to the orginal assumptions. This 
proves the desired result concerning trk. Note that it follows form this result 
that |tr/c| has a strictly positive lower bound on the given interval. 

Estimates for the spacetime will be obtained with the help of a quantity Z. 
This is modelled on the quantity S introduced for similar purposes in [7j which 
in turn was influenced by the arguments of |13j . 

Z = 3i/2 _ Snfip) = 87:p + ^a'^a,, - ^R. (27) 

Here the equality of the second and third expressions follows from the modified 
Hamiltonian constraint (llSp . The time derivative of Z is given by 

— = -2HZ - 2H [(Jijfj'^ + 47r(p + trS*)] < -2HZ. (28) 

Since it has been proved that _ff > it can in particular be concluded that Z 
is non-increasing. Now introduce a quantity l{t) which is positive and satisfies 
l/l = H . This is a replacement for a in the isotropic case. At a big rip singularity 
I tends to infinity. It follows from that 

|(logZ)<-2|(logO. (29) 

Hence d/ dt{\og{ZP)) < 0. As a consequence Z < Ci/~^ for a positive constant 
Ci. It follows that p < (Ci/87r)?"2 ^^lis can be put back into the modified 

Hamiltonian constraint to obtain a useful lower bound for l/l, using the fact 
that in Bianchi types I- VIII the scalar curvature R is non-positive. Compare 
the remarks following equation ()23|) . There results an inequality of the form 
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I > (72?"''"^, where the positive constant C2 depends only on Ci and B. This in 
turn means that must tend to zero within a finite time only depending on 
C2, provided the solution exists that long. Thus it has been shown under the 
given assumptions that if every solution continues to exist as long as / remains 
bounded then every solution will end in a big rip singularity in finite time. The 
existence statement is dependent on the choice of matter model and will be 
addressed later. 

It will now be shown that the geometry isotropizes in the approach to the 
big rip singularity in these models. This means by definition that aija^^ /H"^ 
tends to zero as t tends to the time where / tends to infinity. Since Z is 
bounded by a constant times as t ^ it follows from the definition of 
Z that aij(j^^ — 0{l~'^) in this limit. On the the other hand by the modified 
Hamiltonian constraint it follows that > ^Bl~^'^ and so the quotient of 
interest decays at least as fast as as t — > t*. The quotients p/H"^ and 

R/H^ satisfy similar estimates which says that the density becomes negligi- 
ble compared to the critical density and the spatial scalar curvature becomes 
unimportant in the approach to the big rip singularity. 

It is possible to refine these statements so as to obtain more detailed in- 
formation about the asymptotics of the geometry in the limit t —> t^t. This 
is analogous to the refinement of Wald's theorem proved by Lee [6] . The first 
step is to determine the leading order asymptotics of tr/c. From the modified 
Hamiltonian constraint and the estimates already obtained it follows that 

(trfc)2 = 247r/(p)(l + o(l)). (30) 

Consider now the right hand side of equation (pO| . 

R + {trkf + 47rtr5 - 127rp = {trkf{l + o(l)). (31) 

For a power law Cardassian term pf'{p) = nf{p). In the case of a perfect fiuid 
with equation of state p — {"f — l)p the explicit contribution containing trS 
coming from pS)) adds to AT:pf'{p) to give 4:n-fpf'{p) — 4:TTnjf{p). Hence 

12^(-2/(p) + {p + tYS)f'{p)) = 127r(n7 - 2)/(p) = ^{nj - 2)(trfc)2(l + o(l)). 

(32) 

With the information available concerning Z the estimates previously obtained 
for the terms containing Vij' and k"^iS^m can be improved. Putting all these 
facts together shows that 

|(trfc) = ^(trfc)2(l + o(l)). (33) 

It follows that trfc = ;^(t* - t)^"^(l + o(l)) and that / = lo{t^, - t)^ for a 
constant Iq. Thus the exact leading order asymptotics for trk and I have been 
obtained. 

With trfc under control the leading order asymptotics of the metric can also 
be determined. The starting point is the equation -j^igij) = —^guk'-j. It follows 
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from this that ^{l "^Qij) — —21 "^gua^ j. By using some Unear algebra as in [S] 
it is possible to obtain an estimate of the form 

Wl-^m < \\l-'9{to)\\ + f Ms)\\\\l-^g{s)\\ds. (34) 

J to 

Here to is some fixed time and the norms are the matrix norms of the matrices 
of components of the corresponding geometric objects. The norm of a can be 
estimated in terms of aija^^ for which a decay estimate in the approach to 
t = has already been obtained. There result bounds for the norm of l~^g and 
for the components l~^gij. An analogous argument can be applied to bound 
the components l^g^-' ■ Thus the metric l~^gij is uniformly equivalent to the flat 
metric with components Sij. Using the evolution equation for l~^gij again shows 
that this quantity converges to a limit g*j as i — > i,. Its inverse also converges 
so that g*j is non-degenerate. Finally the leading term in the expansion of 
the metric can be derived. The results which have been proved assuming a 
continuation criterion will now be summed up in a theorem. 
Theorem Consider a spatially homogeneous solution of Bianchi type I- VIII of 
the Cardassian equations pS]) and (fT9)l in the case that f{p) = Bp" with B > 
and n < 0. Suppose that the matter is described by a perfect fluid with equation 
of state p=(7 — 1<7<2. If the solution is expanding at some time then 
the determinant of the metric tends to infinity at a finite time in the future. 
For t 

g,, (t) ^ iU - t)^g:, + o((t, -t)^) (35) 

and trfc — > — . 

71-y 

Note how similar this result is to that in the homogeneous case. As in that 
case more detailed asymptotics can be obtained if desired [T . To complete the 
proof of this theorem a continuation criterion must be proved. Is is contained 
in the following Lemma. 

Lemma Under the assumptions of the theorem if the solution exists on a time 
interval [to,ti) for some finite ti and the determinant of the metric remains 
bounded as f — > ti then the solution can be extended to an interval [to,t2) for 
some t2 > ti. 

Proof A key fact in proving this lemma is the boundedness of Z . The energy 
density satisfies the evolution equation 

^ = (trfc)(p + p)- V,f + a'.kh = ti-kjp + 0{p). (36) 

here it has been used that Z is bounded in order to bound the term containing 
a'j. This implies that 

|(log(p/37))^0(l). (37) 

It can be concluded that since / is bounded on the given interval the same is true 
of Using the fact that Z is bounded again shows that trfc is bounded. It is 
already known that and p arc bounded. The bounds for trfc and aijcr''^ can 
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be used to bound Qij, g^^ and kij by a method which was already used above. It 
follows from the boundedness of trfc that (trfc)~^ is also bounded, as was shown 
following . From this point the rest of the proof proceeds essentially as in 
section 4 of [TU] . The solution remains in a compact subset of the region where 
the coefficients of the evolution equations for the second fundamental form are 
regular as a consequence of the fact that trfc cannot approach zero. 

The analysis of the lemma can be extended to other reasonable matter mod- 
els. This includes mixtures of non-interacting fluids with linear equations of 
state. It also includes the case of collisionless matter satisfying the Vlasov 
equation. This provides a model of the distribution of galaxies which includes 
velocity dispersion. The only new element which comes in is the need for a suit- 
able local existence theorem with continuation criterion. This can be proved in 
the same way for these matter models as in the case of a single fluid already 
discussed. The proof follows the arguments in section 4 of [10] as before and, in 
the case of collisionless matter, those of section 2 of [9]. This implies a partial 
analogue of the theorem for the other matter models. The statements about the 
occurrence of a big rip singularity and the isotropization of the geometry as it 
is approached go through without change. More precise information concerning 
the leading terms in the asymptotics depend on the details of the matter model 
chosen. For it is necessary to know the leading order behaviour of the pressure 
as a function of the energy density. In the case of a mixture of fluids the fluid 
with the smallest value of 7 will dominate. In the case of collisionless matter 
the leading order behaviour will be like that of dust since for collisionless matter 
tiS/ p — > when the scale factor tends to infinity. 

4 Discussion 

In this paper the Cardassian models in cosmology were extended to a large 
class of geometries which are homogeneous but not necessarily isotropic. It was 
shown that for Bianchi types I- VIII, matter described by a perfect fluid with 
linear equation of state and Cardassian terms of the form Bp^ with B > 
and n < the property of isotropic models that all solutions develop a big rip 
singularity extends to all the homogeneous models. Moreover the homogeneous 
models approach isotropic and spatially flat models in a strong sense as the big 
rip singularity is approached. 

At the same time it was pointed out that there is no obvious way of extending 
further to the inhomogeneous case so that it is not really possible to call the 
Cardassian models a theory of gravitation. This statement does depend on 
the point of view that reinterpreting the homogeneous and isotropic models in 
terms of an exotic fluid and passing to the inhomogeneous case in the context 
of the Einstein equations coupled to that fluid is not a satisfactory solution to 
the problem. Arguments supporting this point of view were presented. It is 
difficult to rule out that there is some natural general inhomogeneous version 
of the Cardassian models but the authors have failed to find one. 

At a time when new models designed to address the question of the nature of 
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the cause of accelerated expansion are proliferating, some of them differing from 
general relativity, it would be good to pose the question at an early stage of the 
history of such a model whether there is indeed a general theory behind it. There 
are cases where this is clear, as with the f{R) theories which ar presently very 
popular or other higher order theories such as those with Lagrangian RapR"^ . 
This includes the fact that the field equations for these theories are known to 
be well posed [8] . The example of the Cardassian models shows that things are 
not always so simple. It may be that the range of applicability of a proposed 
model does not go beyond the homogeneous case. 

Big rip singularities can occur for other models which have been introduced 
to explain accelerated cosmological expansion. This applies in particular to gen- 
eral relativity coupled to certain types of exotic matter. It may well be that the 
isotropization phenomenon exhibited in this paper for Cardassian models also 
occurs in the approach to big rip singularities produced by other mechanisms. 
This remains to be investigated. 

5 Appendix: motion in a potential 

Let y be a smooth function defined on an interval (a_,a+) of real numbers 
where a_ and aj^ may be finite or infinite. To exclude trivialities suppose that 
V is not constant. Consider the equation a = V'{a). The energy ^ + V{a) is 
independent of time. Corresponding to initial data at some time for a solution 
of this equation there is a solution on a maximal interval t+) beyond which 
it cannot be extended smoothly. Here t- and t+ may be finite or infinite. If 
one of them is finite then as t tends to this endpoint a{t) must come arbitrarily 
close to a_ or a+. For if not conservation of energy implies that a remains 
bounded. Then a and a remain in a compact subset of the admissible values 
where the coefficients of the equation are regular and by a standard result on 
ordinary differential equations [5J the solution can be extended, contradicting 
the original assumption. 

Let a{t) be a solution with initial data given by a(to) = ao a-nd d(io) = ai- 
Its energy \s E — ^af + V{ao). Let &_ — inf{a : V{a) < E on (6_,ao)} if the 
set over which the infimum is to be taken is non-empty and 6_ = otherwise. 
Correspondingly, let 6+ = sup{a : V{a) < E on (ao,6+)} if the set over which 
the supremum is to be taken is non-empty and = oq otherwise. If 6_ = &+ 
then V'{ao) = and the solution is time independent. Otherwise the solution 
defines an interval [b- , 5+] of positive length. Suppose that > a_ and that 
the solution reaches the point 6_ at a finite time ti > to- Then d{ti) = and 
V{b-) < 0. If V'{b-) were zero the solution could only be the time independent 
solution at 6_. Thus V'{b-) < and a is greater than 6_ at times slightly after 
ti. In summary, if the solution reaches a point 6_ in the interval (a_,a+) 
after finite time it must turn around and start increasing. The same argument 
applies to the evolution backwards in time. An analogous statement can be 
obtained with b- replaced by by a similar argument. A solution behaves in 
a monotone way as long as it is strictly between 6_ and 6+ since a can never 
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vanish there. Thus if in a given solution d(to) > is a will increase for all t > to 
unless b+ e (a_,a+). In the latter case a will change sign and a will start to 
decrease. If a{to) < then a{t) will decrease for all t >to unless b- G (a_, a+). 
Similar statements apply to the past time direction. If both and 6+ are 
in (a_,a+) then there are two times ti and t2 at which the solution reaches 
b-. At these times the initial data are identical and so the solution is periodic, 
oscillating between the two endpoints and Thus the following cases are 
possible. The solution is time independent (type 1). The solution is periodic 
(type 2). The solution is such that a changes sign exactly once (type 3). The 
solution is such that the sign of a is constant (type 4). Which of these types 
occur is easily decided in terms of the position of 6_ and b+. In types 1 and 2 
this already determines the essential features of the dynamics. In types 3 and 
4 it remains to decide the following. How does the solution behave during the 
monotone approach to the point a_ or a+? It can be assumed without loss of 
generality that a{t) > in the time interval of interest and approaches a+. It 
should be remembered that a+ may be finite or infinite and that, independently 
of this, i+ may be finite or infinite. For the case being considered here to occur 
it must be the case that V{a) < E close to a+. By conservation of energy 
d = \/2{E — 2V{ay). Thus what happens depends on the details of the limiting 
behaviour of F as a — > a+ and its relation to E. In the examples of relevance 
in this paper it is always the case that a_ = and a+ = oo. 
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